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qq ■ The chiral phase transition is investigated within the framework of the linear sigma model at finite 

^\ ' temperature. We concentrate on the meson sector of the model and calculate the finite temperature 

^\ . effective potential in the Hartree approximation by using the Cornwall- Jackiw-Tomboulis formalism 

of composite operators. The effective potential is calculated for N = 4 involving the usual sigma 
and three pions and in the large N approximation involving N — 1 pion fields. In the N = 4 case 
we have examined the theory both in the chiral limit and with the presence of a symmetry breaking 
term which generates the pion masses. In both cases we have solved the system of the resulting 
, gap equations for the thermal effective masses of the particles numerically and we have investigated 

the evolution of the effective potential. In the N = 4 case there is indication of a first order phase 
transition and the Goldstone theorem is not satisfied. The situation is different in the general case 
using the large N approximation, the Goldstone theorem is satisfied and the phase transition is of 
the second order. For this analysis we have ignored quantum effects and we used the imaginary time 
formalism for calculations. 
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O '■ I- INTRODUCTION 
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The study of matter at very high temperatures and densities and of the phase transitions which take place between 
JZ* the different phases is very interesting from several points of view and it has been the subject of intense study the 
latest years since it is relevant to particle physics, astrophysics and cosmology. According the standard big bang 
O 1 model it is believed that a series of phases transitions happened at the early stages of the evolution of the universe, 
D | the QCD phase transition being one of them Q-gj- There is present hope that it could also be possible to probe 
• the underlying physics of QCD in the laboratory in experiments involving relativistic heavy ions collisions. These 
\ experiments are planned for the near future and the results could be possible to shed some light in questions like the 
• »h , restoration or not of chiral symmetry, the nature of quark gluon plasma and the physics of neutron stars Q . 
y\ ' Our aim is to study the chiral symmetry of QCD which is spontaneously broken by the small current quark masses. 
A powerful method in approaching questions like the restoration of spontaneously broken symmetries is to construct 
order parameters which characterise the state of symmetry of the system under consideration. These quantities are 
zero in the one phase but not in the other. A classic example of an order parameter is the magnetisation of a 
ferromagnetic substance, it is non-zero below the Curie temperature, but it disappears at temperatures higher than 
that. The system undergoes a transition from an asymmetric, ordered state with non zero magnetisation at low 
temperature to a symmetric disordered state with zero magnetisation at high temperatures well above the Curie 
point. We usually encounter two types of phase transitions. In transitions of the first order the order parameter 
"jumps" discontinuously from its value in the one phase to that in the other (usually zero). In contrast during second 
order transitions the order parameter vanishes continously. 

An important order parameter which is related to the chiral phase transition of QCD is the quark condensate, a 
measure of the density of quark-antiquark pairs that have condensed into the same quantum mechanical state. They 
fill the lowest energy state -the vacuum of QCD- and as a result the chiral symmetry is broken, since there is no 
invariance under chiral transformations We expect that if we raise the temperature the quark condensate will 
disappear and the theory will be chirally symmetric. 

As in many other cases in physics, in order to deal with the chiral phase transition we can use effective models 
to describe the physical situation. In the case of chiral symmetry a model with the correct chiral properties is the 
linear sigma model, a theory of fermions (quarks or nucleons) interacting with mesons |J . This model has been used 
extensively as an effective theory in low energy phenomenology of QCD describing the physics of mesons and it is 
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well suited for a study of the chiral phase transition |?],|| . We review briefly the meson sector of the model and make 
contact with the pion phenomenology in section 3. 

In studies of phase transitions the finite temperature effective potential is an important and popular theoretical 
tool. Early stages of applying similar techniques goes back at seventies when Kirzhnits and Linde jlo| were the 
first proposed that symmetries broken at zero temperature could be restored at finite (high enough) temperatures. 
Subsequent work by Weinberg jnj , Dolan and Jackiw jL2| as well as many others resulted a wide adoption of the 
effective potential as the basic tool in such studies. The finite temperature effective potential V(tfi, T) is defined 
through an effective action T(4>) which is the generating functional of the one particle irreducible graphs and it has 
the meaning of the free energy density of the system under consideration. 

A generalised version is the effective potential V(4>, G) for composite operators introduced by Cornwall, Jackiw and 
Tomboulis (CJT) |Q. According to their formalism a generalised version of the effective action is introduced, which 
in contrary to the usual effective action T(cj>) depends not only on 4>{x) but on G{x,y) as well. These two quantities 
are to be realized as the possible expectation values of a quantum field 4>{x) and the time ordered product of the 
field operator T(j>{x)(j>{y) respectively. In this case the effective action T(<j),G) is the generating functional of the two 
particle irreducible vacuum graphs (a graph is called "two particle irreducible" if it does not become disconnected 
upon opening two lines Jl3|). This formalism was originally written at zero temperature but it has been extended at 
finite temperature by Amelino-Camelia and Pi where it was used for investigations of the effective potential of the 
A</> 4 theory [jl4| and gauge theories 

Physical solutions demand minimization of the effective action with respect to both <fi and G JT^,|IJ]. As a result 
the CJT effective potential should satisfy the stationarity requirments 



dV(<j>,G) 



= 



(1) 



and 



dV(^G) 
dG 



. 



(2) 



Then the conventional effective potential results as V(<j>) = V((f>; Go (</>)) at the solution G(tf>) = Go(</>) of the second 
equation. 

There is an advanta ge i n using the CJT method to calculate the effective potential in Hartree approximation. 
According to reference jL4j , using an ansantz for a "dressed propagator" we need to evaluate only one graph that of 
the "double bubble" in fig. fa, instead of summing the infinite class of "daisy" and "super daisy" graphs, given in 
figs, lb, lc respectively, using the usual tree level propagators. 



(a) 





(c) 



FIG. 1. The double bubble and examples of daisy and superdaisy diagrams 

We demonstrate the advantage of this method in the next section where we calculate the finite temperature effective 
potential for one scalar field with quartic self-interaction. The calculation of the effective potential by using the CJT 
formalism is reviewed in great detail in references (fjjl^l but in order to illustrate the method for the calculation of 
the effective potential for the linear sigma model at finite temperature we will reproduce here the basic steps. 

In our calculations we use the imaginary time formalism, known as Matsubara formalism (f^JT^ [HJ. According 
to this technique we work in Euclidean space-time and use the same Feynman rules as at zero temperature but 
evaluating momentum space integrals we replace integration over the time component k^ with a summation over 
discrete frequencies which means that in the case of bosons ki = 2ninT, n = 0, ±1, ±2, ±3, . . .. This is encoded into 
the following relationship 



2 



where j3 is the inverse temperature, (3 — l/fc^T, and as usual Boltzmann's constant is taken Ub = 1. For the shake of 
simplicity in what follows we have introduced a subscript to denote integration and summation over the Matsubara 
frequency sums. 

The outline of the remaining sections is as follows. In next section 2 we calculate the effective potential for the A(f) A 
theory using the CJT method. In section 3 we apply the same technique to calculate the effective potential for the 
linear sigma model in the Hartree approximation and we solve numerically the system of the resultant gap equations 
both in the chiral limit and with the presence of a linear term which breaks the chiral symmetry of the Lagrangian 
and generates the pion observed masses. We have repeated these steps for the generalised version of the linear sigma 
model with N — 1 pion fields in the large N approximation. Finally, in section 4 we outline our results and conclusions. 



II. THE A$ 4 THEORY 



The Euclidean Lagrange density for a single scalar field with quartic self interaction interaction is given by 



C = ±(d li *)(&**)-±m*&- ^<& 4 . (4) 



where the "mass squared" m 2 is considered as a negative parameter, in order to realise the spontaneous breaking of 
symmetry. By shifting the field as $ — > <f> + </> the "classical potential" takes the form 

u(4>) = \m 2 4> 2 + ^ A . ( 5 ) 

and the interaction Lagrangian which describes the vertices of the shifted theory is given by 

£nt(<r;0 = ~^ 3 -^$ 4 . (6) 
The tree-level propagator which corresponds to the above Lagrangian density is 

V^ifok) = k 2 +m 2 + ^A0 2 . (7) 

According to CJT formalism JT3] the finite temperature effective potential is given by 

V(ct>,G) = U(<t>) + ±jlnG- 1 ( c t>; k) 
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where U(4>) is the "classical potential" given by equation (0) and Vi{4>, G) represents the infinite sum of the two- 
particle irreducible vacuum graphs. We are going to evaluate the effective potential in the Hartree approximation 
which corresponds to that the leading contribution to last term Va {4>, G) of the effective potential comes from the 
"double bubble" diagram given in fig. la [0,0]. Therefore the effective potential results as 



V{<j>, G) = V + ^A0 4 + i J In (T 1 ^; k) 

1 C 1 1 r r 1 2 



+- / [(k 2 +m 2 + -A0 2 )G(^; k) - 1] + -A 



G(^;fc) 



(9) 



Minimizing the effective potential with respect to "dressed propagator" G(cf>; k) we obtain a gap equation 



G _1 (0; fc) = k 2 + m 2 + ^A^ 2 + Ja / G(^; fc) . (10) 

The solution Go(4>;k) of the gap equation is inserted back into the expression for the effective potential to give the 
potential as a function of (f>. According to jl3| using the propagator Go(0; k) for internal lines corresponds to summing 
all daisy and super-daisy diagrams using the usual tree level propagators as in O] . 
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In order to proceed and according to reference fL4j , it is convenient to adopt the following form for the propagator 
G(0;fc) 

where an "effective mass" M = M(</>; k) has been introduced. Then the gap equation for the propagator becomes an 
equation for the effective mass 

where it is obvious that in this approximation the effective mass M is momentum independent. 
In terms of the solution Mq(4>) of the gap equation jl^), the effective potential takes the form 

V(4>, Mo) = \vr?e + A 4 + 1 J ln{k 2 + M 2 } 



(13) 

1(3 ' I" •>'.. ° \J : 

Performing the Matsubara frequency sums as in [T^ ] , the logarithmic integral which appear into the above expression 
of the effective potential divides into a zero temperature part Qo(M) which is divergent and one non zero Qp(M) 
temperature part which is finite and can be written as 

Q(M) = \ f ln(fc 2 + M 2 ) = Q (M) + Qp{M) 

1 J/3 

d 3 k ujy- If d 3 k , , , „ , , 

PY^iR 1 ^ 1 -^ 1 ' (14) 

where lo^ — (k 2 + M 2 ) 1 / 2 . In this above expression and in what follows we omit the subscript of M. Similarly, the 
second integral is divided into two parts as well. One zero temperature part Fq{M) and a finite temperature Fp(M) 
part as 

F{M) = I vTW 2 = Fo{M) + Fp{M) 

d 3 k 1 1 r d?k 1 1 

J2^2^ + PJ (2^^exp(-^ k )-l ' (15) 

The second term vanishes at zero temperature, while the first term survives but it gives rise to divergences which can 
be carried out using appropriate renormalisation prescriptions p^JT^ |. If one is interested in temperature induced 
effects only, as it is our approximation, the divergent integrals can be ignored. In this case, by making a change of 
the integration variables the finite temperature part of F(M) can be written as 

_ t 2 r°° x 2 dx i 

^ >~2^J [x 2 + y 2 } 1 / 2 exp[x 2 + y 2 ] 1 / 2 - 1 ' ( ' 

where we have used a shorthand notation and y = M/T. Similarly the finite temperature part of the logarithmic 
integral becomes 

Q (M) = — dxx 2 \n l-cxpH^ + y 2 ] 1 / 2 ) . (17) 



2tt 2 

Then the finite temperature effective potential, ignoring quantum effects, can be written as 

V(cf>, M) = ^m 2 4> 2 + ^A0 4 + Q P {M) 



\{M 2 -m 2 - \cj?)F {M) + ^[Fp{M)f . (18) 



We can obtain a more compact form if we make use of the gap equation (n3) , 



M) = im 2 2 + 1a^ 4 + Qp{M) - l -\ [Fp(M)f . (19) 
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III. THE EFFECTIVE POTENTIAL OF THE LINEAR SIGMA MODEL 



A. The Linear Sigma Model 

As we mentioned in the introduction, the linear sigma model serves as a good low energy effective theory in order one 
to have some insight into QCD. The model was first introduced in the sixties || as a model for pion-nucleon interactions 
and has attracted much attention recently specially in studies involving the Disoriented Chiral Condensates (DCC's) 
p(]| pl| . The model is very well suited for describing the physics of pions in studies of chiral symmetry. Fermions are 
inserted in the model either as nucleons, if one is to study nucleon interactions or as quarks. The mesonic part of 
the model consists of four scalar fields, one scalar isoscalar field which is called the sigma a field and the usual three 
pion fields n°, tt^ which form a pseudo-scalar isovector. The fields form a four vector (c, 7Ti), i = 1,2, 3 which we 
regard as the chiral field and the model displays an 0(4) symmetry. The a field can be used to represent the quark 
condensate, the order parameter for the chiral phase transition, since they exhibit the same behaviour under chiral 
transformations The pions are very light particles and can be considered approximately as massless Goldstone 
bosons. We use the model as an effective theory for QCD ignoring the fermion sector at the moment and concentrating 
on the meson sector. 

The generalised version of the meson sector of the linear sigma model called the O(N) or vector model and is based 
on a set of N real scalar fields. The O(N) model Lagrangian can be written as 

£ = i(^*) 2 -im 2 * 2 -^A* 4 -ea, (20) 

and in the absence of the last term it remains invariant under O(N) symmetry transformations for any N x N 
orthogonal matrix. In order our notation to be consistent with applications on pion phenomenology we can identify 
the $i with the a field and the remaining N ~ 1 components as the pion fields, that is €> = (cr, 7ri, . . . , ttn-i)- The 
last term, ea into the above expression has been introduced in order to generate the observed masses of the pions. 
The contact with phenomenology is obtained by considering the the case N = 4. Then the Lagrangian of the model 
is given by 

C = \(daf + \{d^f - imV - - A (o * + ^ f £a , (21) 

where e = /^m 2 and f n — 93 MeV is the pion decay constant. At zero temperature and in order to be consistent 
with the pion observed mass of m T ss 138 MeV and the usually adopted sigma mass m a « 600 MeV we choose the 
coupling constant A in our model to be 

3(to 2 — ml) , . 

A = V a *' . (22) 

The negative mass parameter m 2 has been introduced in order to obtain spontaneous breaking of symmetry and its 
value is chosen to be 

- m 2 = (to 2 - 3m 2 )/2 > . (23) 

As we have referred earlier, our aim is to use the linear sigma model as a model to study the chiral phase transition, 
so we need to calculate the effective potential. For this calculation we adopt the CJT formalism which has been 
presented earlier for the A</> 4 theory and calculate the effective potential in the Hartree and large N approximations. 



B. Hartree approximation in the chiral limit e = 

In order to deal with the exact chiral limit first, the starting point is the Lagrangian @ and we ignore the 
symmetry breaking term at the moment. By shifting the sigma field as a — > a + <j) the "classical potential" results as 

U&) = \m><? + A 4 (24) 

and the interaction Lagrangian which describes the vertices of the new theory takes the form 

dnM 4) = —a 4 - ^ - A a V - ^ - , (25) 
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plus terms linear in the a field and constants which we omit for simplicity. In our approximation we do not consider 
interactions given by the last two terms in the Lagrangian. 

The tree level sigma and pion propagators corresponding to the above Lagrangian are 



V- 1 ((f>;k) = k 2 +m 2 + ^X<l) 2 . 



(26) 



(27) 



We evaluate the effective potential in the Hartree approximation which means that we only need to calculate the 
"double bubble" diagrams as in X(j> theory. In the linear sigma model the corresponding effective potential at finite 
temperature can be written as 



V(<f>,G) = U{4>) 



lnG-^fc) 



+ - 



1 



[V-\frk)G a {k)-l\ 



V-\<k k)G^k) 



-V 2 {<i>,G a ,G*) , 



(28) 



where the first term U(<fi) is the "classical potential" and the last term 1^(0, G CT , G„) originates from the sum of 
"double bubble" diagrams |Ti}| . There are four types of double bubbles as we show in fig. 2 and these contribute the 
following terms in the potential 



Vv 



, G a , Gjr ) — 3— 



GAfrk) 



4 



G*{4>;k) 



4 



G^-k) 



GAfrk) 



(29) 




3 9 6 6 

FIG. 2. The double bubble graphs which contribute to the effective potential for the 0(4) linear sigma model in the Hatree 
approximation. The numbers show the weight of each type of bubble in the expression for the effective potential. 



Minimizing the effective potential with respect to the "dressed" propagators, we get the following set of nonlinear 
gap equations 

G- 1 (q>-k)=V- 1 (frk) + ^ [ G,(0;fc) + ^ / G CT (0;fc) 



G- 1 ^; k) = V- 1 ^; k) + ^ / G w (& k) + ~ / G CT (0; k) . (30) 



6 J ' 6 

In order to proceed we can use the same ansatz for the dressed propagators as in the one field case JL4| 

G7^=fc 2 +M ff 2 A , (31) 
and using the equations (§fj), ©, @ and (|l|) we end up with the following system for the thermal effective masses 

M 2 a =m 2 + \M 2 + \F{M a ) + ~F(M n ) 

M 2 =m 2 + -\4> 2 + —F(M a ) + ^F(M n ) . (32) 
6 6 6 
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In these last two equations we have used a shorthand notation and F(M) is given by 

As in Xcf) 4 theory the thermal effective masses are independent of momentum and functions of the order parameter <fi 
and the temperature T. 

By using these two equations the effective potential at finite temperature then can be written as 
V(cj>, M) = \m 2 4 2 + ^A0 4 + ~ J ln(fc 2 + M 2 ) + ln(fc 2 + M 2 ) 

\{Ml - m 2 - l^ 2 )F(M a ) - \{Ml -m 2 - \m 2 )F{NU) 

A c; \ \ 

+ -lF(M a )} 2 + -[F(M.)] 2 + —F(M a )F(M n ) . (34) 

Minimizing the effective potential with respect to the "dressed" propagators we have found the set of nonlinear gap 
equations for the effective particles' masses given by equation (p2|). In addition, by minimizing the potential with 
respect to the order parameter we obtain one more equation 

= m 2 + \\(f> 2 + —F{M (T ) + -F(M W ) . (35) 

In order to study the evolution of the potential as a function of temperature, we perform the Matsubara frequency 
sums jijj as in the one field case. There are some troubles concerning the renormalisation of the model [ p^ , ^6|j2^j27| , |28t . 
At the level of our approximation we ignore quantum effects for the moment and keep only the finite temperature 
part of the integrals. Using a compact notation the finite temperature effective potential can be written in the form 

V(<f>, M) = im 2 2 + ^A0 4 + Q p (M) + 3Qp(M v ) 

A ^ A A 

- -lF (M a )} 2 - —[F/3(M n )] 2 - — Fp [M a ) Fp (M n ) , (36) 



where in this last step we have used the gap equations (|32|). The expessions for Fp(M)aiid Qp{M) are given by the 
equations (16) and (|l7j). 

In order to calculate the effective masses as functions of temperature we need to solve the system of equations ( |3^ ) 
and (^||). We first observe that if = 0, which happens in the high temperature phase, the two equations become 
degenerate, the particles have the same mass and we have to solve only one equation 

M 2 =m 2 + XF p {M) , (37) 

where obviously as in the expression of the effective potential we keep only the finite temperature part of the integral. 
This last equation can be used to define a "transition temperature" T cl . This temperature is defined as the temperature 
where both particles become massless. Recall now that Fp(M) is given by 

T 2 f°° x 2 dx 1 

F ^ M ^ = 2^J [x 2 + (M/T) 2 ] 1 / 2 exp[x 2 + (M/T) 2 } 1 / 2 - 1 ' (38 ^ 

where we show explicitly the dependence of Fp(M) to mass and temperature. When the mass of the particles vanishes 
the integral reduces to the well known 

. . f°° xdx it 2 
^=1 ^ (39) 

and the temperature T c \ is found to be 

T cl =V2(-^) 1/2 . (40) 

But in defining our model parameters we have chosen that at zero temperature (f> 2 — f 2 — — 6m 2 / A, where f v = 
93 MeV is the pion decay constant, so we find that T c \ = V^f^ » 131.5 MeV. 
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In the low temperature phase we can eliminate <fi and end up with the following nonlinear system 

Ml = -2m 2 - XF {M a ) - XF (M W ) 
Ml = ~F p (M a ) + ^F I3 (M 7T ) , 

which we solve numerically and the solution is presented in fig. 3a. 



(41) 
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FIG. 3. (a) Solution of the system of gap equations in the chiral limit. The sigma and pion effective masses are given as 
functions of temperature, (b) Evolution of the order parameter as a function of temperature. 

As shown in fig. 3a, the temperature T cl which is calculated numerically was found to be in excellent agreement with 
the value obtained by using the limit of the high temperature equations with degenerate masses. 

At this point we can observe that there is an indication of a first order phase transition, because combining the 
first of the two equations in (|32|) with equation (|3^) we find that 



Mi 



(42) 



This last equation shows of course that the order parameter varies with temperature proportionally to the sigma 
mass. The temperature dependence of cj) = 4>{T) is calculated by using the sigma mass as it was found by solving the 
system in eqn (fi"l|). This is given in fig. 3b where it is obvious that this approximation predicts a first order phase 
transition. This last observation coincides with the qualitative picture given by Baym and Grinstein in their early 
paper p5| where their "modified Hartree approximation" predicts a first order phase transition as well. However, in 
contrast to our approximation into their analysis they included quantum effects as well. Signals of a first order phase 
transition have also been reported in recent analyses by Randrup |^6| and by Roh and Matsui p7| . 

In order to get more insight into the nature of the phase transition and verify that the transition is of the first order 
we can calculate the effective potential V(<f>, T) as a function of the temperature and the order parameter. So, we first 
solve numerically the system of equations (|3^) and (^) (where of course we keep only the finite temrerature part of the 
integrals) and calculate the effective masses of the particles as functions the order parameter and the temperature. 
Finally the effective potential is calculated numerically by using these masses. The evolution of the potential for 
several temperatures is given in fig. 4. The shape of the potential confirms that a first order phase transition 
takes place, since it exhibits two degenerate minima at a temperature T c sa 182 MeV which is usually defined as 
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the transition temperature. The second minimum of the potential at <f> ^ disappears at a temperature T C 2 ~ 
187 MeV. The temperatures T c \ and T C 2 are called (in condensed matter terminology) the lower and upper spinodal 
points respectively. Between these temperatures metastable states exist and the system can exhibit supercooling 
or superheating. For T c \ < T < T c the metastable states are centered around the origin since for T c < T < T C 2 
the metastable states occur for <f> ^ 0. When the system reaches T cl or T C 2 the curvature of the potential at the 
metastable minima vanishes. A discussion about first order phase transitions and more details about how these 
transitions procceed can be found in r cfs ||-§]. 
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FIG. 4. Evolution of the effective potential V(4>, T) as a function of the order parameter <j) for several temperatures in steps 
of 2 MeV. The two minima appear as degenerate at T c ps 182 MeV. 

C. Hartree approximation in the broken symmetry case e 7^ 

When e 7^ the term linear in the sigma field into the Lagrangian generates the pion observed masses. This term 
is G independent and so minimization of the potential with respect to "dressed" propagators will give us the same 
set of gap equations for the effective masses as before. However, minimizing the potential with respect to <fi we get 
the following equation 



±\4? + ^F p {M a ) + -Ffi(M v ) 



'-£ = 



(43) 



In order to proceed we need to solve the nonlinear system of three equations (|3 
T = the above equation becomes 



and (p 



We first observe that at 



Ml 



I A 2 = - = ml 
6 V 6 



(44) 



where is the tree level pion mass. Then for <f> — f n we recover the relation between the pion mass at zero 
temperature and the symmetry breaking factor e: e = f^m^ where / w is the pion decay constant. We solved the 
system of eqns ( |32| ) and ( [43| ) numerically and the solution is presented in fig. 5a. At low temperatures the pions 
appear with the observed masses but their mass increases with temperature since the sigma mass decreases. At high 
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temperatures (higher than ~ 300 MeV) due to interactions in the thermal bath all particles appear to have the same 
effective mass. 

The appearance of the symmetry breaking term into the Lagrangian modifies the evolution of the order parameter 
(j) as well. As it is obvious in fig. 5b, as the temperature increases the order parameter decreases and at very high 
temperatures vanishes smoothly. But in this case the change is not a phase transition any more. We rather encounter a 
smooth crossover from a low temperature phase when the particles appear with different masses to a high temperature 
phase where the thermal contribution to the effective masses makes them degenerate. 




FIG. 5. (a) Solution of the system of gap equations in the case when t / 0. At low temperatures the pions appear with the 
observed masses, (b) Evolution of the order parameter as a function of temperature. 



D. Large N approximation in the chiral limit e = 

The large N approximation of the linear sigma model has been studied recently by Amelino-Camelia ]28fl and our 
expressions are very similar to the ones obtained there, since the same method is used in both cases. However, in our 
approach we do not consider the renormalization of the model because in our approximation we take into account 
only finite temperature effects. Our analysis is in a sense complementary to that in ps|| since we solve the system of 
gap equations and consider the effects of the symmetry breaking term (the last term in the Lagarangian (EG)) which 
is omitted in reference |2g] . 

The Lagrangian of the linear sigma model when we consider a large number (N — 1) pion fields is given in equation 
( pp| ) and shifting the sigma field as a — > a + <fi, the tree level propagators are 

o \ 

D- 1 (cj ) :k) = k 2 + m 2 + — 4? , (45) 



D-\d > ;k) = e+m 2 + — 4> 2 . (46) 



Then the effective potential at finite temperature will appear as 
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;fc)G CT (0;fc)-l] + 



N — 1 



[p- 1 (0;fc)G 7r (0;fc)-(7V-l)] 







where the last term originates from the double bubble diagrams and its contribution reads as 



V 2 (4>, G a , Gtt) = 3 



A 

6iV 



~] 2 



G a {k) 







X(N 2 - 1) 
6iV 



I 2 



X(N-l) 
67V 



G CT (0;fc) 



GAfr k) 



(47) 



(48) 



The weight factors appearing in the above expression can be understood in a similar way as in the 0(4) case the 
only difference being the N — 1 pion fields. Of course it is easy to see that we recover the previous case by simply 
substituting N = 4. 

As in the case of Xcj) 4 and the 0(4) model we minimise the effective potential with respect to the dressed propagators 
and we get a set of gap equations. By using the same form for the dressed propagators as before, we end up with the 
following set of nonlinear gap equations for the thermal effective particle masses 



2X 



2A 
~N 



Fp{M a 



2\(N - 1) 
3N 



Ml = m 2 + + + 



(49) 



where we only keep the finite temperature part of the integrals. As it is easy to observe, for N = 4 we obtain identical 
expressions for the system of gap equations as in the case of the 0(4) model. Then the effective potential will appear 
in the form 



y(^,M) = imV + ^A0 4 + i 



A 



2N 



X(N 2 - 1) 



\n(k 2 + Ml) 
2A 



Fp(M„)f 



X(N — 1) 
3N 



F p {M a )F f) {M v ) 



(50) 



In the large N approximation, which means that we ignore terms of 0(1/N), the system of the two equations (|4£ 
reduces to 



Ml 
Ml 



2X i7 2X , , 
— <* 2 + ~^Fp{M n ) 



N 
3N 



2X 



Fp{M„) 



(51) 



We have retained the terms quadratic in <f> since <j> depends on N as <p 2 = 3Nm 2 /2X and so these terms are of O(l). 
In order to solve this system and be in "some contact" with phenomenology in the chiral limit, we set now N = 4 so 
at zero temperature the pions are massles and sigma has a mass Ml = —2m 2 . Now our system is written as 



1 2X 
Ml = m 2 + -X4> 2 + —F p (M 7r ) 

1 2A 
Ml = m 2 + -A0 2 + —F P (M«) 



(52) 



and in order to solve it we proceed in a similar way as in Hartree approximation. 

At very high temperatures the potential has only one minimum that at <f> — and in this case the two equations 
become degenerate 



Mi = M. 



J-2 



M 



2A 



Fp(M) 



(53) 



This last equation actually defines the critical temperature. Fp{M) is given by the same expression as in the 0(4) 
case. The mass of the particles vanishes at the critical temperature, so we can use the result of eqn. (^) to find that 
the critical temperature is at 
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6m 



2 \ V2 



V3f n w 161MeV 



(54) 



Before proceeding to examine the low temperature phase we should note an observation which actually marks 
the significant difference between the Hartree approximation in the N — 4 case and the large N approximation. 
Minimizing the potential with respect to <p gives 



dV{<l>, M) _ dV 



= 



which in the large N approximation becomes 

dV(4>,M) _ 
10 ~ 



m 2 + ^ + ^F (M„] 



= 



Combining this last equation with the second of equations (52) above we observe that 



dV{(f>, M) 



<pMl = 



(55) 



(56) 



(57) 



Therefore the large TV approximation implies that the pions should be massless in the low temperature phase in 
accordance with the Goldstone theorem. 

This observation is reflected in the solution of the system of the gap equations as it is shown in fig. 6a. The pions 
at low temperatures appear as massless, but at high temperatures the thermal contribution to the effective masses 
make them degenerate with the sigma. The order parameter vanishes continuously in this case as it is shown in fig. 
6b and this corresponds to a second order phase transition. 
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FIG. 6. (a) Solution of the system of gap equations in the large N approximation in the chiral limit. At low temperatures 
the pions appear as massless. (b) Evolution of the order parameter with temperature in the large N approximation. 



E. Large N approximation in the broken symmetry case £ / 

As already mentioned for the 0(4) case, the symmetry breaking term ea has been introduced into the Lagrangian 
in order to generate the observed masses of the pions. The same can be done for the O(N) model the only difference 
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being the N — 1 pion fields. Inserting this term into the expression for the effective potential and differentiating with 
respect to <j> we obtain, as in the 0(4) case, one more equation. In the large N approximation this is written as 



m 2 + iA0 2 + ^(F,(M.) 



'-£ = 



(58) 



We have solved this last system of equations (|32|) and ( |58| ) numerically and the solution is given in fig. 7. As in 
the N = 4 case there is no longer any phase transition. We encounter again the crossover phenomenon between the 
low and high temperature phases, the difference being now that the change of the order parameter (fig. 7b) in the 
transition region is much more smooth in contrast the more "sharp" behaviour seen in the N = 4 case in fig. 5b. 




FIG. 7. (a) Solution of the system of gap equations in the large N approximation in the case of broken chiral symmetry 
e / 0. At low temperatures the pions appear with the observed masses, (b) Evolution of the order parameter with temperature. 



IV. CONCLUSIONS 



We have studied the chiral phase transition using the linear sigma model. In order to get some insight into how 
the phase transition could proceed, we have calculated the finite temperature effective potential of this model in 
the Hartree and large N approximations using the CJT formalism of composite operators. The method proved to 
be very handy since we actually need to calculate only one diagram. In both cases we have solved numerically the 
system of gap equations and found the evolution with temperature of the effective thermal masses. In the Hartree 
approximation we find a first order phase transition, but in contrast to that, the large N approximation predicts a 
second order phase transition. This last observation seems to be in agreement with different approaches to the chiral 
phase transition based on the argument that the linear sigma model belongs in the same universality class as other 
models which are known to exhibit second order phase transitions |t|,|2^| . However in the large N approximation the 
sigma contribution is ignored and this of course introduces errors when we calculate the critical temperature. In the 
case N = 4 which is closer to phenomenology we could probably obtain a better approximation if we had considered 
the effects of interactions given by the last two terms in the Lagrangian ( p5| ) . We are planning some investigations in 
order to include the effects of these terms in the calculation of the effective potential. 

When we include the symmetry breaking term ea which generates the pion observed masses, both in Hartree 
and large N approximations we found that there is no longer any phase transition. We rather observe a crossover 
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phenomenon where the change of the order parameter in the Hartree case occurs more rapidly in contrast to the more 
smooth behaviour exhibited by the large N approximation. Our observation confirms results reported recently by 
Chiku and Hatsuda |3^] using a different approach. In their analysis they also report indication of first order phase 
transition in the chiral limit. 

Of course as we have already pointed out, the linear sigma model is only an approximation to the real problem 
which is QCD, but the study of the chiral phase transition in the framework of this model could be a helpful guide to 
how one could tackle the original problem and get some insight in the physics invloved. We have used the imaginary 
time formalism which is adequate for studies at thermal equilibrium but if one is interested in studies of the dynamics 
of the phase transition, the real time formalism seems to be more convenient |3|,[l8|,[l9| . We are currently making some 
preliminary investigations in this direction. 
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